Nuclear symmetry energy and core-crust transition in neutron stars: a critical study 
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The slope of the nuclear symmetry energy at saturation density L is pointed out as a crucial 
quantity to determine the mass and width of neutron-star crusts. This letter clarifies the relation 
between L and the core-crust transition. We confirm that the transition density is soundly correlated 
with L despite differences between models, and we propose a clear understanding of this correlation 
based on a generalised liquid drop model. Using a large number of nuclear models, we evaluate 
the dispersion affecting the correlation between the transition pressure Pt and L. From a detailed 
analysis it is shown that this correlation is weak due to a cancellation between different terms. The 
correlation between the isovector coefficients -fsTsym and L plays a crucial role in this discussion. 
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Stimulated by the development of exotic nuelear 
physics, the efforts to determine the nuclear equation of 
state (EOS) have focused in the last few years on the den- 
sity dependence of the symmetry energy S{p) In 
particular, the symmetry-energy slope at saturation den- 
sity, represented by the quantity L, has raised a great 
deal of interest [l|-ly|: while the different nuclear models 
widely disagree on the value of this basic quantity, in- 
creasing experimental data are expected to bring 
more and more stringent constraints, leading to a radi- 
cal progress in our knowledge of the EOS of neutron-rich 
matter. Several experimental constraints on the slope 
L have already been proposed in the last decade. Such 
studies include information obtained from the mass for- 
mula 0, isospin diffusion Q, experimental double neu- 
tron to proton ratio [91, and isoscaling parameters in 
heavy ion collisions |10|, pygmy dipolc resonances 
giant dipole resonances |12| . neutron-skin thickness |5[- 
This impacts strongly on the physics of compact stars. 
In this letter, we will discuss the link between L and 
the transition from the liquid core to the solid crust of a 
neutron star. It has been claimed that a precise determi- 
nation of L would give a tight indication of the density 
Pt and pressure Pt at the transition point , and conse- 
quently the mass and extension of the crust which play 
a crucial role in the interpretation of pulsar observations 
[l3| . However, the role of L in the determination of the 
core-crust transition needs to be checked against model 
dependence and clarified, as mentioned in Ref. In 
the present work, we use a variety of nuclear models to 
address this issue. We verify and explain the strong cor- 
relation between L and pt- However, we show that when 
independent models are considered there is no real corre- 
lation between L and the pressure at the transition point. 
This behaviour results from a competition between op- 
posite effects which destroy the correlation. This seri- 
ous limitation has to be taken into consideration when 
drawing astrophysical consequences from the experimen- 
tal determination of L. 

Catalyzed matter in compact stars satisfies the /3- 
equilibrium condition which favors very neutron-rich 
matter: the proton fraction is reduced to a few percent 
in the region of the core-crust transition. Compact-star 



structure crucially depends on the symmetry energy, for 
a wide density range. The density dependence of the 
symmetry energy, S{p), is deduced from the energy den- 
sity functional obtained in the framework of mean field 
nuclear models. Besides, it can be expressed as a develop- 
ment around the saturation density pq, whose coefficients 
correspond to the isovector parameters of a generalised 
liquid-drop model (GLDM): 
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where x = {p — po)/{3po). Here and in the sequel, 
the index "IV" ("IS") attributed to the coefficients of 
the GLDM stands for "isovector" ("isoscalar"). The 
first coefficients have received traditional denominations: 
civ,o = J = S{po), civ,i = L, civ,2 = Ksyra, etc. In the 
framework of the parabolic approximation, the energy 
per particle for asymmetric matter is given by E(p, y) = 
E{p,0) + S{p)y'^, where y = (/?„ — Pp)/p- For conve- 
nience, we will use in the following either the isospin- 
asymmetry y or the proton fraction Yp — {1 — y)/2. 
This approximation allows to emphasise the role of the 
GLDM coefficients, so we will use it to analyze our re- 
sults, although the calculations have been performed us- 
ing the complete density functional of each model. In 
the parabolic-GLDM framework, the energy per particle 
reads: 



Eip,y) 



n>0 



ClS,i 



ClV,ny^) 



(2) 



In the isoscalar channel, we have cis,o = Eq = E{po), 
cis,i = 0, cis,2 = Koo, etc. 

We will show results obtained from a set of non- 
rclativistic and relativistic effective interactions, together 
with results from a microscopic Brueckncr-Hartree-Fock 
(BHF) calculation using the interaction Avl8 [l^ with 
Urbana three-body forces [l^. As non-relativistic effec- 
tive models, we take Skyrme-type interactions from dif- 
ferent groups (SV, sen, RATP, SkMP, Gs, Rs, SkI2, 
SkI3, SkI4, SkI5, SkI6, SlylO, Sly230a, Sly230b, Sly4, 
SkO, NRAPR, LNS, BSkl4, BSkl6^Skl7); the respec- 
tive references can be found in [l7H2G |. Besides, we con- 
sider two different types of relativistic effective nuclear 
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models: (i) non-linear Walecka models with constant 
couplings (NL3 HiLtMI [13, GMl, GM3 H, FSU, 
NLw/O [23|, NL/5(5 [23); (ii) hadronic models with den- 
sity dependent coupling constants (TW DD-MEl, 
DD-ME2 113, DDH,5 0). Let us remark that the EOS 
features present more variability within the relativistic 
models than within the Skyrme ones [l7l | . 

The inner crust of a neutron star is usually modelised 
as a lattice of very neutron-rich nuclei, immersed in a 
gas of electrons and dripped neutrons. As the density in- 
creases, the difference between the nuclei and surround- 
ing neutron gas decreases, until the stellar matter be- 
comes homogeneous: this is the transition to the liquid 
core, when the pasta phase dissolves [Hlflll^- In order 
to determine the transition point, one should in principle 
compare the free energy of homogeneous matter to that 
of any inhomogeneous configuration, which constitutes 
the pasta phase [2§|. It has been shown in [s^l that the 
composition of the crust, namely the size of the pasta 
structures, is sensitive to the very low density neutron 
matter EOS [3l| while almost not affecting the crust-core 
transition. Therefore, we expect that our conclusions will 
not depend on the very low density behavior of the neu- 
tron matter EOS. In [U, [T^ it has been verified that the 
transition density obtained from the pasta phase calcula- 
tion can be very well approximated by the entrance into 
the dynamic spinodal region, under the constraint of (3 
equilibrium. The spinodal is the density region where 
the homogeneous matter is unstable against density fluc- 
tuations, due to the nuclear liquid-gas phase transition 
affecting the bulk EOS. In the case of finite size density 
fluctuations, the Coulomb and surface terms reduce the 
instability: the dynamic spinodal region is then smaller 
than the thermodynamic one obtained when only the 
bulk term is considered. The difference between the two 
regions is in principle model dependent, via the nuclear 
surface term. However, we found that this effect is too 
small to play a role in the present discussion. Thus, for 
simplicity, wc will focus on the transition density pt, pro- 
ton fraction Yp.t and pressure Pt taken at the crossing 
point between the /3 equilibrium EOS of stellar matter 
and the thermodynamic spinodal, keeping in mind that 
they represent shifted values of the actual density, proton 
fraction and pressure at the core-crust transition. 

It has been noticed in previous works that the transi- 
tion density decreases as L increases [H 0, IB] ; this cor- 
relation has been verified with many different models, 
see Fig. [1^). We have found that this behaviour can 
be understood through the energy-density curvature of 
pure neutron matter (NM), denoted Cnm- More specif- 
ically, we have considered the curvature Cnm taken at 
the density of the upper spinodal border in symmetric 
matter, p^. This quantity, denoted Cnm,s = CnmIPs) ~ 
^ ^rfpiT"^ {Ps), where E^m is the energy per particle in 
neutron matter, is indeed correlated with the transition 
density pt, as shown in Fig. [2^). This result allows a qual- 
itative interpretation. The spinodal region corresponds 
to the region of {p, Yp) where the energy density has a 
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FIG. 1: (Colour online) Correlation between L and a) the 
transition density pt, b) the transition proton fraction l^,t, 
c) ifsym. The full (empty) symbols are for Skyrme forces 
(relativistic models) and the asterisk for BHF. 



negative curvature. For very asymmetric matter such 
as /3-equilibrium matter, Cnm,s gives a good indication 
to localise the position of the spinodal border: the larger 
Cnm,s is, the farther should be the spinodal contour from 
the point (p = Ps,Yp = 0); and this corresponds to a 
lower Pt- Besides, Cnm.s is strongly correlated with L: 
see FiglSb)- This relation appears clearly when Cnm,s 
is expressed in the parabolic approximation, in terms of 
the isovector coefficients civ.n: 
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Note that the isoscalar terms of this expansion are ex- 
actly zero at p = Ps by definition of ps- Since x is neg- 
ative, the influence of the higher order terms n > 2 is 
weakened. Furthermore, for all the models considered, 
we have ps — (2/3)po: this makes the contribution of the 
term n = 2 in Eq. ^ close to zero. As a result, Cnm,s de- 
pends very weakly on i^sym and is essentially determined 
by L, as can be verified on the figure. In summary, the 
correlation observed between L and pt can be understood 
as the consequence of the link existing between L, Cnm,s 
and Pt- 

The proton fraction at the transition, Yp^t, is also ex- 
pected to decrease with increasing L. Indeed, a smaller 
symmetry energy corresponds to a lower proton fraction 
for /? equilibrium matter. Assuming a consensual value 
of J (about 32 MeV) , a larger L means a smaller symme- 
try energy at subsaturation densities and, consequently, 
a smaller l^,t, as is shown in Fig. [Ud). The dispersion 
of data in this figure reflects the model dependence of J . 
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FIG. 2: (Colour online) Correlation of the energy-density cur- 
vature of neutron matter Cnm.s with a) the transition density 
pt, and b) the symmetry energy slope L. The full (empty) 
symbols are for Skyrme forces (relativistic models) and the 
asterisk for BHF. 



We also show in Fig. [T};) the correlation between ifsym 
and L, which will be useful for the following analysis. 

Let us now turn to the study of the transition pressure 
Pt . In order to investigate what evolution of Pt we should 
expect for increasing L, we express the pressure in the 
parabolic-GLDM framework: 



P{p,y) 
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We can first notice that, for fixed density and asymmetry, 
the pressure should increase with L. For instance, within 
the parabolic GLDM, the pressure of neutron matter at 
saturation is just Pnyiipo) = Pipo, 1) = Lpo/3. At sub- 
saturation density and in neutron matter, there is still a 
correlation between L and Pnm that we show in Fig. [5Ji) , 
for the density p=0.08 fm^'^. An equivalent result was 
obtained in [l[ at p^O.l fm~'^. In asymmetric matter 
corresponding to the transition position (pt, yt), this sim- 
ple relation is spoilt by two effects: the dependence of 
the transition position (pt, yt) on L and the dependence 
of the isovcctor coefficients civ.?i appearing in the pres- 
sure ([U on L. We indeed remind the correlation between 
(pt, j/t = 1 — '^Yp^t) and L shown in Figs. [T^)-b): pt de- 
creases with increasing L while yt increases with increas- 
ing L. From Eq. ([3]), these additional correlations might 
change the evolution of the pressure Pt with respect to 
L. In Fig. we represent Pt versus L calculated 
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FIG. 3: (Colour online) The pressure versus L a) for neutron 
matter at p = 0.08 fm~^; b) for /3-equiIibrium stellar matter 
at the transition point. The full (empty) symbols are for 
Skyrme forces (relativistic models) and the asterisk for BHF. 



consistently for each of the models considered. If we con- 
sider only the sub-group formed by the Skyrme models, 
it could be noticed a very slight decreasing correlation 
of Pt for L >60 MeV. However, considering all the mod- 
els in Fig. [U there is a large dispersion for the values of 
Pt, and it is not possible to define a general trend. We 
can notice that the dispersion of Pt {L) corresponds to an 
amplification of the dispersion already present in PnmI^) 
at p~0.08 fm^'^. The dispersion in Fig. means that 
the pressure is quite sensitive to the GLDM coefficients 
others than L: thus, it is difficult to extract a clear re- 
lation between Pt and L which would be satisfied by all 
the models. 

Let us now analyze the result shown in Fig. In 
the parabolic-GLDM, the pressure (jSj depends clearly 
on various terms which are correlated to L. In the fol- 
lowing, we analyze the contribution of these correlations 
to the global relation between Pt and L. The variation of 
the transition pressure Pt with L can be decomposed into 
two different contributions: the first one is related to the 
variation of Pt with respect to the transition position, 
{pt,yt)i and the second one is related to the variation of 
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FIG. 4: (Colour online) Expected variation of Pt with L, com- 
paring the diflterent contributions: variation due to the posi- 
tion, position 1 (down-triangles) and position 2 (up-triangles) , 
and due to civ.n coefficients, coef 1 (diamonds) and coef 2 
(squares). The resulting total pressure variation is very close 
to zero (circles). The full (empty) symbols are for Skyrme 
forces (relativistic models). 



Pt with respect to the isovector coefficients in Eq. 
civ,n- We now introduce our notations: the variation of 
the transition position induces two contributions: 
(hereafter named position 1), and (named posi- 

tion 2). The variation with respect to the isovector coeffi- 
cients in the GLDM induce the terms ^£iv^ (named 



coef n). All the partial derivatives are taken at the tran- 
sition position (pt,yt) obtained with the complete func- 
tional for each model. The variation of pt, yt and Kgym 
with respect to L is extracted from linear fits in Figs.[TJi)- 
c), giving the values 5pt/5L = (-3.84 ± 0.24) x 10"^ 
MeV-i fm-3 and 8yt/5L = (6.08 ±0.82) x lO^^ MeV-^, 
and SKsyn-i/SL = 3.07 ± 0.33. Notice that in doing so, 
we smooth out the dispersion that still exist between the 
different models (see Figs.[T^)-c)). However, anticipating 
the results, the dispersion is the largest for the variable 
Yp^t which has a very weak contribution to the global vari- 
ation of Pt. We consider the contributions of the terms 
coef n up to n = 2. We have verified that no general 
correlation appears between L and civ,n for n > 2; fur- 
thermore, due to the fast convergence of the expansion 
([3]) at p = Pt, the contribution for n > 2 is small. 

We show in Fig. 2] the contributions of the terms po- 
sition 1, position 2, coef 1 and coef 2 to the variation of 
Pt- The contribution coef 1, namely ^ = p^y'^/{3po), is 
clearly large and positive, as expected. This term alone 
would predict an increase of the transition pressure with 
L, which is not observed. However, this term is balanced 
by the sum of the two terms position 1 and coef 2, which 
are negative and of similar magnitudes (down-triangles 
and squares in Fig. 2]) . The term position 2 brings only 
a negligible contribution. The sum of the different con- 
tributions, represented by the circles in Fig. 21 is thus 



very close to zero; furthermore, its sign depends on the 
model. We conclude that there is no clear correlation be- 
tween L and Pt, due to the non- trivial balance between 
opposite contributions which differ among models. This 
can explain why opposite behaviors have been obtained 
for Pt{L) in the literature [3, Hi]. 

It is important to notice the role of the L — Kgym cor- 
relation, leading to the contribution coef 2: 
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If this term had not been considered, the prediction that 
Pt increases with L would have persisted, since the term 
position 1 is not sufficiently large to overcome the term 
coef 1. In other words, the decrease of pt is not sufficient 
to predict a decrease of Pt if we consider only the leading 
term of the GLDM development in Eq. ([3]). The correla- 
tion L — A'syiii is therefore crucial in this discussion. 

In the present letter we have explained why there is a 
good correlation between the crust-core transition den- 
sity and the symmetry energy slope, L; and we predict 
that this behaviour should not depend on the relation 
between L and i^sym- On the contrary, no correlation of 
the transition pressure with L was obtained. We have 
highlightened the competing contributions to the varia- 
tion of Pt with L: this explains that the effect of L on the 
transition pressure is not strong enough to overcome the 
dispersion due to the interaction details of independent 
models. This means that an experimental determination 
of L alone will not be sufficient for a good estimation 
of the crust mass and moment of inertia of a compact 
star. In fact, the range of variation of Pt obtained in 
the present letter with a large set of nuclear models lies 
within the interval indicated in [l3|, 0.20 < Pt < 0.65 
MeV/fm"^, but completely out of the interval obtained in 
0, where a correlation between Pt and L was supposed 
and the range of L values was determided from isospin 
diffusion. The large dispersion of the predicted transition 
pressure obtained when independent models are consid- 
ered needs to be reduced. A more accurate knowledge of 
the isoscalar EOS could improve the situation. Exper- 
imental constraints given directly at subsaturation den- 
sity for the isovector EOS would also help. The striking 
correlation between ifsym and L appears to be an impor- 
tant feature, which should be further investigated. 
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